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Abstract 

The neutron-rich n Li halo nucleus is unique among nuclei with known separation energies by its 
ability to emit a proton and a neutron in a f3 decay process. The branching ratio towards this rare 
decay mode is evaluated within a three-body model for the initial bound state and with Coulomb 
three-body final scattering states. The branching ratio should be comprised between two extreme 
cases, i.e. a lower bound 6 x 10~ 12 obtained with a pure Coulomb wave and an upper bound 
5 x 10 -10 obtained with a plane wave. A simple model with modified Coulomb waves provides 
plausible values between between 0.8 x 10~ 10 and 2.2 x 10~ 10 with most probable total energies of 
the proton and neutron between 0.15 and 0.3 MeV. 
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I. INTRODUCTION 



Among their remarkable properties, nuclei with a neutron halo display unusual /3-decay 
channels. There is indeed the possibility that the decay occurs in the halo, releasing the 
halo nucleons. This process has been observed in the (3 delayed deuteron decay of 6 He and 
It is however severely limited by the energy conservation condition 



ii 



Li 



S 2n < 5( 2 H) + (m„ -m p - m e )c 2 « 3.007 MeV (1) 



where S2 n is the two-neutron separation energy of the halo nucleus, 5( 2 H) is the binding 
energy of the deuteron, and m n , m p , and m e are the neutron, proton, and electron masses, 
respectively. Only a few nuclei have low enough separation energies to allow this decay: 
6 He, 8 He, n Li, 14 Be, 17 B, 27 F, . . . 

Another /3 decay process is even more severely restricted, i.e. a decay of a halo neutron 
releasing a free neutron and a free proton. The condition is 

S 2n < (m n -mp- m e )c 2 « 0.782 MeV. (2) 

Among nuclei with known two-neutron separation energy, the unique nucleus where this 
decay is allowed is 11 Li, 

n Li ^ 9 Li + n+p + e~ + i> e (3) 

with the separation energy [7] 

S 2n = 378 ± 5 keV. (4) 

This process should be observable if the branching ratio is large enough. However, the 
small energy available for the decay indicates that the phase space is much smaller than for 
the deuteron emission. How rare is this decay is the main question raised in the present 
exploratory study. 



The n Li nucleus is described in a 9 Li+n+n three-body model [8| as in our previous 
studies of the deuteron delayed emission 0, [loj]. The 9 Li+n+p final state is in the three- 
body continuum of n Be. The calculation of wave functions in this continuum is much more 



complicated than in the three-body continuum of 6 He [111, [12| . The construction of three- 
body scattering states for 9 Li+n+n would already be more difficult than for a+n+n because 
of the poor knowledge of the 9 Li+n interaction. The study of the 9 Li+n+p continuum is 



worse for several reasons, (i) The halo nucleons are not identical and the wave functions 
have about twice as much components at the same level of truncation in an expansion 
in hyperspherical harmonics, (ii) The presence of a Coulomb interaction between the 9 Li 
core and the proton requires a more complicated treatment than in the neutral case, (iii) 
The structure of the continuum wave functions is more complicated since one can expect 
a larger number of bound states to which they must be orthogonal. For these reasons, 
the technique that we have developed [l2| can not provide a converged calculation with our 
present computer capabilities. Since an evaluation of the branching ratio would be necessary 
to guide future experiments, we shall simplify the study by describing the continuum with 
three-body Coulomb waves. This approximation should be accurate enough to estimate the 
order of magnitude of the branching ratio and the shape of the energy distribution. 

In section |TT1 we present general formulas for the decay probability per time unit for the 
ft delayed np emission. In section IHIl we evaluate the branching ratio and discuss its origin. 
Concluding remarks are presented in section IIVI 



II. DECAY PROBABILITY FOR (3 DELAYED np EMISSION 
A. General expression of decay probability 

In this section, we establish the general expression for the f3 decay probability distribution 
for a three-body final state in the continuum. The initial nucleus with mass number A is 
described as a three-body bound state of a core and two nucleons. This state with angular 
momentum Jj, projection Mi, and parity 7Tj is expressed in hyperspherical coordinates. The 
spin, isospin and parity of the core are neglected. Three-body scattering states are discussed 



in Ref. 



121 ] and used in Ref. [13]. 



Let us follow the notation in Ref. 13] and denote the three particles as 1, 2, and c (for the 
core). Let k\ 2 be the relative wave vector between particles 1 and 2 and fe c (i2) be the relative 
wave vector between the center of mass of those particles and the core. When the spin of the 
core c is neglected, the outgoing scattering states can be denoted as , , where 

M\ and M 2 are the projections of the spins of particles 1 and 2. These states are assumed 
to be normalized with respect to 5(ki 2 — k' 12 )5(k c (i 2 ) — k' c ^)5 M\_m'J> m 2 m^- 



The distribution of decay probability per time unit can be written as 
dW 1 m e c 2 



GifiQ - e)^— £ £ (W + a 2 e i m ^i 2 I , (5) 

% Mi M X M 2 \ l* 



dki 2 dk c (i2) 2tt 3 H 13 

where Gp ~ 2.996 x 10~ 12 is the dimensionless /3-decay constant, A ~ —1.268 is the ratio 
of the axial-vector to vector coupling constants, and E is the total energy of the nuclear 
fragments. The Fermi integral f(Q — E) depends on the kinetic energy Q — E available for 
the electron and antineutrino with 

Q = (m„ - m p - m e )c 2 - S 2n . (6) 

The Fermi and Gamow- Teller matrix elements are respectively given by 

i=i 

and 



M GT ,(E) = 2<*£ fc J £ V^|*^-> (8) 

where Sj and tj are the spin and isospin of particle j, and /x = —1, 0, +1 labels the tensor 
components of the spin. 

If one integrates expression (JSJ) over all directions, the distribution of probability as a 
function of the total center-of-mass energy E < Q of the three particles is given by 

" W - 1 ^G lm -E)(?m + X* d -^). (9) 



dE 2tt 3 h ""' J 

The Fermi and Gamow- Teller strengths appearing in this expression are given for o = F or 
GT by 

dB(a)_ i r JU r JU j„ m 2 m 



do) 



Mi MiM 2 M 



dE 2.7, + lJ "J ^ y 2/i 12 2/i c(12) 

where \i\ 2 is the reduced mass of particles 1 and 2, and // c (i2) is the reduced mass of the core 
c and the system 1+2. 

The total transition probability per time unit W is obtained by integrating (Q from zero 
to Q. The branching ratio can than be derived as 

TZ = Wt 1/2 /\n2, (11) 

where t±/ 2 ~ 8.75 ms is the half life of n Li. 



B. Bound-state and scattering three-body wave functions 



In hyperspherical coordinates, the three-body wave function of a bound state is defined 

as 

* JMw (p, n 6p ) = p~ 5/2 Yl 4k(p) y J ,K^ P i (12) 

■yK 

where p is the hyperradius, fl^n represents the five hyperangles, and y^xi^p) l $ a hyper- 
spherical harmonics (see Ref. [8] for definitions and notations). The symbol K corresponds 
to the hypermomentum quantum number and 7 is a shorthand notation for l x l y LS, where L 
is the total orbital momentum, S is the total spin, and l x and l y are the orbital momenta for 
the relative motions corresponding to the Jacobi coordinates between particles 1 and 2 and 
between the core and the center of mass of 1+2, respectively. The parity of these relative 
motions is given by 7T = (— Y) lx+l v = (— 1) K , which implies that K must be either even or 
odd. The hyperradial wave functions are expanded as 

N 

xiK(p) = J2 c %Mp) (13) 

in terms of iV Lagrange functions fa (see Ref. [8] for definitions). Since the hyperspherical 
harmonics and the Lagrange functions are orthonormal, the wave functions are normed if 

N 

EEte) 2 - 1 - (14) 

-yK i=l 

In the present approximation of n Li with a + core, particles 1 and 2 are neutrons. The 

effective angular momentum and parity are J 77 = + . The isospin is T = 1 for the halo 

neutrons. Since they are identical, antisymmetry imposes (— l) lx = (— l) s . 

The final states of the decay are three-body scattering states. It is convenient to replace 

the projections Mi and M 2 by the total spin 5* of nucleons 1 and 2 and its projection v. 

With a zero-spin core, 5 is the channel spin. The ingoing scattering states read [131 ] 

3/4 

\ 11 \ 11 it a ?i it t 71 it\ ~\ iL,,,M— v * ir, 

lj XUJl'yUJ -'-'UJ -t^LU 

x E(- 1 )'"^(^)X 7 J F( 7 „^)(P) ! (15) 

lK 1 — 1 

where A c = A — 2 is the core mass number. This formula differs from Ref. 13) because of a 



/ A \ 3/4 

*fci3*c(u)s« = pirv 5 ' 2 f^-J E E (L u sM- v »\jM)y£X:z(n t 

V C/ JM l xu l„ u L u K w 



different normalization. The normalization for the hyperradial partial waves is 12] 



X^kJp) "> i K " +l (2*/kf 2 [H- K+2 {kp)5^J KKu - Ujl JuKu H+ K+2 (kp)] . (16) 



p— >oo 



In this expression, the wave number k is given by y/2rriNE/h 2 , where is the nucleon 
mass, and Ufy Ku is an element of the infinite-dimensional collision matrix. The subscript 
to refers to the entrance channel. Let us recall here that, in a three-body scattering state, 
there is in principle an infinity of degenerate entrance channels. 
For charged systems, one has 

H yK+2( x ) = g k+I (VjK,x) ±iF K+ i(rj yK ,x), (17) 

where G K+3 / 2 an d F K+3 / 2 are the irregular and regular Coulomb functions, respectively jl^ |. 
The Sommerfeld parameters t\^k are given by 

where K is a diagonal element of the effective-charge matrix and depends thus on the 
channel. One usually neglects non-diagonal terms of this matrix 15) . 

In the neutral case rj lK = 0, Eq. ffTTj) reduces to an expression independent of 7, 

H± K+2 (x) = ±i [—) [J K+2 (x) ± iY K+2 (x)} , (19) 

where J n (x) and Y n (x) are Bessel functions of first and second kind, respectively. 



C. Reduced transition probabilities 

For 11 Li, with the spin of the core neglected, we assume Jj = Mi = 0. The final state 
is a three-body 9 Li+n+p scattering state (fT5l) . Selection rules restrict this state to its + 
and 1 + components for the Fermi and Gamow- Teller transitions, respectively. In the present 
approximation, the properties of the final state only depend on the total spin S and isospin 
T of the two nucleons. For the nucleons in the 9 Li+n+p continuum of n Be, the isospin is 
given by (— iy^ +s+T = — l. For S = 0, l x even corresponds to T = 1 and l x odd to T = 0. 
For S — 1, l x even corresponds to T = and l x odd to T = 1. The number of channels is 
thus about the double of the number of channels in the 9 Li+n+n continuum of n Li. 

The sum over Mi and M 2 in Eq. (JS]) can be replaced by a sum over the channel spin 
equal to S and its projection v. If one replaces the wave vectors fc 12 and fe c (i2) by their hy- 
perspherical counterparts k and f2 5 fc [3] , the reduced transition probabilities can be written 



as 

dB{a) 



dE 



l^ 2 (^) 3 E/^5 fe |M CT ,| 2 , (20) 



where \i = for F and \i = — 1, 0, 1 for GT. 

After integration over f2s p , the matrix elements can be written as 



^ x uj ^Vuj 



where the spherical harmonics depend on the hyperangles characterizing the wave vectors, 
i.e. they depend on the directions of emission of the core and nucleons, and on the repartition 
of the total energy E between these particles [131 ] . The expressions l/ a ^i Vu L u SK u { (J ) are one ~ 
dimensional integrals over the hyperradius p. After integration over and summation 
over the projections p and u, the reduced transition probabilities simplify as 
dB(a) _2J + l (2m N \ 3 ^ Jw a 

dE ~ (27T)6 \ V 2^ ^Jy^SK^) ■ W 

l xtjj l yuj -'-'to >J J^uj 

Let us list the possible cases. For the Fermi operator, the scattering-state partial wave 
has J = and tt = +1. One obtains for 5 = 0, 

X°i x i x om{i Xui iy w ook u ) (p) Xi x i x ow (p) dp, (23) 

i x even i\ 

and for 5 = 1, 



Xi x i x iiK{i Xu iy u iiK w ){p)Xu x iiK(P)dp- (24) 

i x oaa i\ 

For the Gamow- Teller operator, the scattering-state partial wave has J = 1 and 7r = +1. 
One obtains for 5 = 0, 

IixJ yu wKw(GT) = J- ^2 I x}j x iQK(i x j y ^iMu)(p)Xu x iiK{p)dp, (25) 

V '/. odd K J " 

and for 5=1, 

POO 

IixJvuLvIkS^) = ~ ^2 ^2 Xi x i x oiK(i Xu i Vu L u iK ai )(p)Xi x i x ooK(p)dp 

L x even K 

~\ <? 2-^ 2-^ I Xl x l x nK{l Xw l Vu L ul lK u ) (p)xli x nAp)dp. (26) 
V 6 l x odd K J ° 

Because of the properties of Lagrange functions, the integrals are simply given by 

POO 

/ X J ,K^ K JP)X° + Mdp « ^(/ i A,) 1/2 c^xS (7 ^ u) (to i ), (27) 
Jo 

where Xi and A, are the zeros and weights of the Gauss quadrature associated with the 
Lagrange functions and h is a scaling factor providing mesh points = hxi adapted to the 
extension of the physical system. 



D. Coulomb-wave approximation 



As mentioned in the introduction, we shall use a simpler approximation based on three- 
body Coulomb functions. In the pure Coulomb case, the scattering partial waves are ap- 
proximated as 



X%^K U ){P) = 2i K (27r/k) 5 / 2 F K+3/2 (r )lK ,kp)5 11 J KKu . 
With this approximation, the reduced transitions probabilities become 



dB(F) _ Am N x - * - 
dE ~ irkh 2 ^ 2s 



lm,., Ku. 



where 7^ represents here l x J Xui SS with (— l) s = (— , and 

dB (GT) _ dB(F) 
dE ~ dE 

The F and GT reduced transition probabilities are then proportional. 



(28) 



(29) 



(30) 



III. RESULTS AND DISCUSSION 
A. Q value and Fermi integral 

With the separation energy (TjJ of 11 Li, the Q value for the (3 delayed np emission is quite 
small, 

Q w 0.404 MeV. (31) 

Moreover, the wave number is also small, 

k < 0.14 far 1 . (32) 

This will affect the behavior of wave functions at small distances. 

In Fig. [U the Fermi integrals f(Q — E) for the emission of the different hydrogen isotopes 
are compared. The emitted electron being much faster than the heavy particles, the charge 
Z = 4 is used in the electron attraction by the final nuclear system. The Q values are 2.63 and 
4.82 MeV for 2 H and 3 H, respectively. Both processes have been observed experimentally. 



In spite of a much larger Fermi integral, the branching ratio for tritons [ly, [17| is not larger 
than for deuterons (4 , 5 ] . The emission of deuterons can be fairly well described in a model 

8 



where the 9 Li+d resonance observed in the model of Ref. JldH is shifted to about 0.8 MeV 
and an absorption towards other open channels is included |l8j. To our knowledge, no model 
description of the (5 delayed triton emission is available. The difficulty comes from the fact 
that this decay can not be described in a three-body model. 
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FIG. 1. Fermi integral f(Q — E) as a function of the total energy E of the emitted nuclear fragments 
for the hydrogen isotopes 1 H, 2 H, and 3 H. 



The Fermi integral for /3 delayed proton-neutron decay is much smaller than the other 
ones because of the limited phase space. The branching ratio can thus be expected to be 
much smaller than for the other (3 delayed emissions. 

B. Bound-state and Coulomb wave functions 

The 11 Li ground state is obtained as a bound state in a 9 Li+n+n model. The Minnesota 
force is used as nucleon-nucleon interaction jl9j. The 9 Li-neutron interaction is the P2 



interaction of Ref. 



201 ]. downscaled by a factor 0.97 to approximate the experimental binding 



energy. The s and p3/2 forbidden states are eliminated by supersymmetric transformations 
21] . The sum over partial waves in the wave function (1121) is restricted to K < K max = 20. 



This wave function involves 66 components. The hyperradial functions (fl3|) are expanded 
over a Lagrange-Laguerre basis with integrals calculated with the corresponding Gauss- 



9 



Laguerre quadrature, as explained in Ref. The calculation is performed with N = 40 
basis functions and mesh points and the mesh is scaled with a factor h = 0.4 (see Ref. 8] 
for definitions). The resulting energy is —0.391 MeV, close to the experimental value. We 
use the experimental Q value f[3"Tj) in the calculation of Fermi integrals. 

Because of the low values (132 p of the wave number k, the hyperradial scattering wave 
functions should be small at distances where the integrals (p4|) -( l27|) are significant. They 
become smaller and smaller with increasing hypermomentum K. For this reason, the sum 
in Eqs. (122|) or (|29|) is strongly dominated by K = 0. The K = component of the ground 
state thus plays a crucial role. One should however not expect to use a low value of K max 
because the convergence of this K = component is slow [8]. 

For the three-body scattering states, we use approximations based on Coulomb waves. 
Let us first evaluate the effective charges entering the Sommerfeld parameter. Because of 
the K = dominance in expression (129|) . we can restrict ourselves to this value and average 
the Coulomb potential over the K = hyperspherical harmonics. The Coulomb potential 
is simply 

Vc = r^- r (33) 
\r c - r x \ 

where subscript 1 corresponds here to the proton. Using the hyperradius p and the hyper- 
angle a js|, the K = average can be written as 



P 



-00 
00 



3e 2 



1 00 



p cos a 

3e 2 16 r /2 . 2 

'Pel / sin a cos ada, 

P 7T J 



(34) 



where p c \ = A c / (A c + 1) is the reduced mass of the core and the proton and 7 = represents 

lx = l y = L = S = 0. Hence, the effective charge reads 

48 



y0+ 
^00,00 



4.83. (35) 

7TV 10 

To simplify a calculation dominated by K = 0, we shall use this effective value for all partial 



waves. 



C. Distribution of decay probability per time unit 

Various approximations of the distribution of decay probability per time unit for the 
f3 delayed np decay of n Li are displayed in Fig. [2j With the effective charge fl35l) . one 
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0.1 0.2 0.3 0.4 

E (MeV) 

FIG. 2. Distribution of decay probability per time unit for the /3 delayed np decay of n Li: plane 
wave (upper dashed curve), Coulomb wave with effective charge ([35]) (lower dashed curve), and 
shifted Coulomb waves with a = 10 (lower full curve) and 15 fm (upper full curve). 

obtains the lower dashed curve giving the total probability W = 5.1 x 10~ 10 s _1 and thus 
the branching ratio 1Z = 6.5 x 10~ 12 . These results can be contrasted with a plane-wave 
calculation (r] 00 = 0) which leads to the upper dashed curve giving W = 3.8 x 10~ 8 s _1 
and 1Z = 4.8 x 10~ 10 . The Coulomb-wave calculation is pessimistic because it neglects 
an enhanced probability of presence of the emitted nucleons at short distances due to the 
attractive nuclear interaction. The plane-wave calculation overestimates the probability of 
presence of the emitted proton at short distances because of the missing Coulomb repulsion 
by the nucleus. Both calculations neglect a possible absorption towards other open channels 
affecting the final wave function. However, it is difficult to figure out whether one of these 
cases is a better approximation. Hence we turn to a slightly different approach. 

For a better but still simple approximation based on Coulomb functions, we have con- 
sidered the eigenstates of the 9 Li+n+p system. With K max = 20, its wave functions involve 
121 components. The conditions of the calculation are the same as for 9 Li+n+n except for 
the additional Coulomb interaction (133]) between 9 Li and p and a reduced symmetry. The 
9 Li+p relative motion only requires the elimination of an s forbidden state. 

11 
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FIG. 3. K = components of the lowest positive-energy pseudostate near 0.4 MeV for K m&x = 12 
(dashed line), 16 (dotted line), and 20 (dash-dotted line) normalized to the K = Coulomb wave 
with ^oo^oo = 4-83 (right full line); same Coulomb wave shifted by 15 fm (left full line). 

We obtain four bound states at -12.027, -3.944, -0.876, and -0.786 MeV with respect 
to the 9 Li+n+p threshold. Even the lowest bound state is far above the experimental ground- 
state energy —20.14 MeV. The state at —0.876 is the isobaric analog of the n Li ground state. 
The lowest positive-energy state is located at 0.379 MeV. It must not be considered as a 
resonance but rather as a pseudostate, a bound-state approximation of a scattering state at 
this energy. Its wave function will be useful to construct a better exploratory approximation. 

The K = components of the lowest positive-energy state located near 0.4 MeV obtained 
with K m3X = 12, 16, and 20 are displayed in Fig. [3J The energies do not vary much with 
Kmax. but the wave function is not yet converged. The amplitudes of the three curves are 
normalized to the K — Coulomb wave (128]) corresponding to the charge (1351) (right full line). 
One observes a significant shift between the pseudostate and the Coulomb wave. As a simple 
qualitative approximation, we shift the Coulomb wave by 15 fm towards shorter distances 
(left full line). The resulting curve simulates the general behavior of the pseudostate. A 
shift by 10 fm would also be plausible. 

We thus simulate the K — component of the scattering state with the shifted Coulomb 



12 



functions 



xiK {M (p) = 2t K (2ir/k) 5 / 2 F K+3/2 [ VyK , k(p + a)]5„J KKLj (36) 

with a = 10 and 15 fm. These functions do not vanish at the origin but this drawback has 
little influence, i.e. a smaller influence than other approximations. The results are displayed 
as full curves in Fig. [2j the lower curve corresponds to a = 10 fm and the upper curve 
corresponds to a = 15 fm. Their maximum is slightly shifted towards higher energies. The 
most probable total energies E are located between 0.15 and 0.3 MeV and the most probable 
total energies of the proton and neutron should approximately lie in the same interval since 
the 9 Li core is heavier. This approximation corresponds to 0.6 x 10~ 8 < W < 1.8 x 10~ 8 s _1 
and 0.8 x lO" 10 <TZ<2.2x lO" 10 . 



IV. CONCLUSION 

In this paper, we evaluate the order of magnitude of the branching ratio for the delayed 
np emission by n Li, a very exotic decay process, unique among nuclei with known two- 
neutron separation energies. We have established the theoretical formulas for the Fermi and 
Gamow- Teller transitions leading to three-body final states. 

An accurate model calculation is made very difficult by the need of three-body scattering 
states involving three different particles, two of them charged, at very low energies and by 
our lack of knowledge of physical properties of this three-body continuum and of absorption 
effects in the final three-body channel. To circumvent these difficulties in an exploratory 
calculation, we have made several simplifying approximations. Simple models of the final 
state involving a plane wave and a pure Coulomb wave provide likely upper and lower 
bounds of the branching ratio, respectively. We think that more reasonable estimates of 
the branching ratio and of the energy distribution of the decays are obtained with shifted 
three-body Coulomb functions. 

The obtained branching ratio should be comprised between 6 x 10~ 12 and 5 x 10 -10 with 
more plausible values between 0.8 x 10~ 10 and 2.2 x 10~ 10 . The most probable total energies 
of the proton and neutron should lie between 0.15 and 0.3 MeV. In any case, the branching 
ratio is much smaller than for the deuteron and triton channels, i.e. (1.3 ± 0.13) x 10~ 4 [6j 



and (0.93 ± 0.08) x 10" 4 |lj, respectively. It is even much smaller than for the hindered 

13 



deuteron decay of 6 He, (2.6 ±1.3) x 10~ 6 [3[. The main cause of this smallness is the small Q 
value of the process which leads to a limited phase space. The observation of this (3 delayed 
decay mode, if it is possible, will thus require high radioactive beam intensities and long 
measurement times to reach a significant enough number of n Li decays. 

If this unique decay process is studied experimentally, a better model calculation will 
become necessary, with a full calculation of the three-body 9 Li+n+p continuum wave func- 
tions, using the formalism developed in Sec. [Ill This study should be performed with 9 Li+n 
and 9 Li+p optical potentials in order to take absorption effects into account. 
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